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Abstract 

The exact numerical diagonalization of the Hamiltonian of a 2D circular 
quantum dot is performed for 2, 3, and 4 electrons. The results are compared 
with those of the perturbation theory. Our numerical results agree reasonably 
well for small values of the dimensionles coupling constant A = ^ where o 
is the dot radius and as is the effective Bohr radius. Exact diagonalization 
results are compared with the classical predictions, and they are found to be 
almost coincident for large A values. 
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I. INTRODUCTION 



There is considerable interest in the basic science and technological applications of quan- 
tum dots [1-5]. The electron motion in quantum dots are confined to a region with dimensions 
comparable to the de Broglie wavelength of the particle. The result is the quantization of 
energy. However, since the quantization in the vertical direction is much stronger than in 
the planar directions, a quantum dot can well be treated as a 2D disc of finite radius. 

Although quantum dots and natural atoms contain comparable number of electrons, 
their electronic properties can have gross differences due to the nature of the confinement 
potential [6]. In quantum dots confinement potential is experimentally controllable [4], and 
usually parabolic and much shallower than the - potential of the natural atoms. One recalls, 
however, that parabolic structures are not the only ones because hard wall type confinement 
can also be obtained by etching techniques [7-10]. Electron electron correlations in quantum 
dots have been worked out for parabolic and hard wall confinements in the classical limit 
and finite temperature in [11], and for hard wall confinement at zero temperature in [12]. 

In this work we perform an exact numerical analysis of a few electrons confined in a 
circular 2D quantum dot. We believe that such analyses of 2D structures can shed light on 
the nature of the electron correlations in quantum dots. 

The organization of the paper is as follows. In Sec. 2 we formulate the second quantized 
Hamiltonian and specify its symmetries. In Sec. 3 we perform a numerical diagonalization 
of the Hamiltonian obtained its energy levels as a function of the dot dimension. We also 
present a comparative analysis of the quantum mechanical and classical energy levels. 



II. MODEL AND THE METHOD OF CALCULATION 

We consider a many electron system in a two-dimensional circular quantum dot with 
radius a. The system is described by the Hamiltonian 

where m* is the electron effective mass and e is the dielectric constant of the medium. The 
second term in the paranthesis, t/(x), is the hard- wall confinement potential which van- 
ishes for |x| < a, and infinite for |x| > a. For convenience, we introduce the dimensionless 
Hamiltonian 7i via the definition 

where A = a/ as is the dimensionless coupling constant. The dimensionless Hamiltonian H is 
given by 

N 1 X 1 

i=i z z ijtj l x * x jI 
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One notes that energy is measured in units of ^ , 2 which depends on A, and length is 

TYV CL *\ 

measured in units of a, the dot radius. 

In the second quantized language, the dimensionless Hamiltonian takes the form 

T~C = ^2 £-K<^K a K + n Vk,L,M,N 0'KClM a LO'N , (4) 

K 1 K,L,M,N 

where £k is the single particle energy level, and 

V k ,l,m,n = J /d 2 5d 2 xVx(x)^(x)^(x-xO^(x / )^(x') (5) 

is the matrix element of the Coulombic interaction term in (3). Here v?a(x) (A = K, L, M, N) 
is the single particle eigenstate of the free Hamiltonian, and A is a collective index designating 
the radial (tia), angular momentum (ttia) and spin (a a) quantum numbers. Solution of the 
eigenvalue problem for the free Hamiltonian yields 

V?a(x) = 0„ AimA (x)x CTA , (6) 

where Xa is the spin wavefunction, and nm (x) is the normalized orbital eigenfunction given 
by 

1 1 

n , m (x) = — — ^e imf? J H (A; niH |x|) , (7) 

V 71 " \ J\m\+l{k n ,\m\)\ 

with the eigenvalues £k = \Ki K ,\m K \ which is independent of both spin and the sign of m. 
Here k n ^\ m \ are the zeroes of the Bessel function (J\ m \(k nt \ m \) = 0). 

In determining the spectrum of the total Hamiltonian in (3), it would be convenient 
to find a representation for the Coulombic potential V(5i — x') = r^z^ri in terms of Bessel 
functions, as the single particle eigenstates in (7) do depend only on the Bessel functions. 
Accordingly, we use the following decomposition for l/|x — x'| 

00 POO 

£ / dke^ -^J m (kp)J m (kp')e- kiz >- z <\ (8) 



x — X' 



where we let z > — z< — > as the vertical dimension of the dot is vanishingly small. 

By inserting the equations (6), (7) and (8) into equation (5), we obtain the final form of 
the spin independent potential matrix element Vk,l,m,n as follows 



yn K ,riL,nM,nN 
v m K ,m L ,mM,m N 



J\m K \+l(kn K ,\m K \)\ \J\m L \+l(k nL ,\m L \)\ I J\m M \+l \kn M ,\m M \ 

dk I dppJ\ mK \{k nKi \ mK \p)J\ mL \{k nL> \ mL \p) 



|^|mjv|+l(^n J v,|m JV |)| 

Jm K -m L (kp) / dppJ\ mM \(k nMt \ mM \p)J\ mN \(k nNj \ mN \p)J mN - mM (kp). (9) 

J 

Integration over the azimuthal angle 9 in (5) drops the sum over m in (8), and the 
difference between the orbital angular momentum quantum numbers becomes the order 
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of the Bessel function. In (9) J mK -m L and J mN - mM arise from such angular integrations, 
and unlike the other Bessel functions in (9), they have negative or positive order index. As 
J-\ mK - mL \ = {-l) mK ~ mL J\ mK -m L \ and J_j mjv _ mM | = {-l) mN ~ mM J\ mN _ mM \, these two Bessel 
functions create sign differences among different potential matrix elements. 

From the symmetries of the problem it is seen that the z component of the total spin S z = 
J2iLi s i and the z component of the orbital angular momentum M = X^Li m i are conserved 
quantum numbers. The diagonalization of the Hamiltonian is done in the single particle basis 
by taking the constraints coming from the conservation of z component of total spin and 
orbital angular momentum.Hence, dimension of the total Hamiltonian matrix depends on 
S z and M. In the next section we shall present the result of the exact diagonalization of the 
many body Hamiltonian matrix for 2, 3 and 4 electrons. 

III. NUMERICAL RESULTS AND DISCUSSION 

In this section we present the A dependence of the ground state energies for 2, 3 and 4 
electrons. For each electron number we consider certain values of S z and M.In the figures 
below we plot the ground state energy of the dimensionless Hamiltonian H in (3) as a 
function of A. 

The dimensionless coupling constant A measures the dot dimension in units of Bohr 
radius. In plotting the figures we shall vary A from to 3. For < A < 1 one can analyze this 
problem using perturbation theory as was already done in [13]. However when A exceeds unity 
perturbation theory methods do not work, and one has to resort to some other technique, 
such as numerical diagonalization. Hence, ground state energies will be obtained by the exact 
diagonalization of the dimensionless many body Hamiltonian. 

In Fig.l depicted is the A dependence of the ground state energy for 2 electrons. Here solid 
curve is for 251 (M = 0, S z = 0), dashed curve for 252 (M = 1, S z = 0), and short-dashed 
curve for 2T1 (M = 1, S z = 1). 

Fig. 2 displays the A dependence of the ground state energy for 3 electrons. In this figure 
curves correspond to the ground state energies of 3-D1 (\M\ — 1, S z — 1/2) (solid curve), 
3L>3 (|M| = 0, S z = 1/2) (short-dashed curve), 3D2 (\M\ = 2, S z = 1/2) (dashed curve), 
and 3Q1 (|M| = 0, S z = 3/2) (dotted curve). 

In Fig. 3 we present the case of 4 electrons. The solid curve is for 451 (M = 2, S z — 0), 
dashed curve for 452 (M = 0, S z = 0), short-dashed curve for 453 (M = 1,S X = 0), dotted 
curve for 4T1 (M = 0,S Z = 1), and dot-dashed curve for 4T2 (M = 1, S z = l). 

As is seen in these three figures, our exact diagonalization results are comparable with 
those ones which are obtained from perturbation theory [13]. However, there are some differ- 
ences between our numerical results and perturbative results. For instance, 252 state in Fig.l 
and 3.D3 state in Fig.2 are energetically more favoured than 2T1 and 3Q1 states, respectively. 
Also in Fig. 3, 452 and 453 states are more stable than 4T1 and 4T2 respectively. 

It might be interesting to compare the exact diagonalization results with those of the 
classical considerations. Modelling the quantum dot by an isolated conducting disc with the 
capacitance C = (2ea)/n one gets the classical interaction energy Efj^ ss - = j, in units of 

t~i 2L j~; where N is the number of electrons. In Fig.4 we present the A dependence of the 

quantum mechanical interaction energies for 251 (solid curve), 3D1 (dashed curve) and 452 
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(short dashed curve) in units of m ^ 2L f-- The dotted curve in this figure corresponds to the 
classical interaction energy. In fact, we see that exact diagonalization results always remain 
below the classical result [13], and approach to it for large A values. One further observes that 
for small A values, quantum mechanical and classical results behave differently, as opposed 
to the case of large A. This can be attributed to the effect of electron-electron correlations, 
which implies that the ground state energy cannot be written as the sum of the contributions 
of kinetic energy and interaction potential [12]. 

Figure 5 shows the dependence of the dimensionless ground state energy per particle, 
So/N, on the particle number N for different A values. For small A values, So/N deviates 
from a linear behaviour due to the dominance of the free Hamiltonian.On the other hand, 
when the interaction potential becomes dominant with increasing A value, So/N changes 
almost linearly with N, as is observed in the figure. 

In conclusion, we have worked out the exact diagonalization of the total Hamiltonian for 
2, 3 and 4 electrons in a 2D circular quantum dot as a function of the dot dimension. Our 
results approximately agree with those of the perturbation theory for small radii, but dis- 
aggrement sets in for larger radii. We have also compared exact diagonalization results with 
the classical predictions, and found that they almost coincide for large A values. 
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IV. FIGURE CAPTIONS 



Figure 1: The ground state energy levels for 2 electrons 
Figure 2: The ground state energy levels for 3 electrons 
Figure 3: The ground state energy levels for 4 electrons 

Figure 4: Comparison of the classical interaction energy with the lowest three quantum 
mechanical interaction energies 

Figure 5: Dimensionless ground state energies per particle Sq/N versus the particle number 
iV for A = 5(0), A = 10(+), A = 50(D), A = 100(x) 
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Fig.l 
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Fig.2 
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Fig.3 
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